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Abstract 

Dynamics of relativistic outflows along the rotation axis of a Kerr 
black hole is investigated using a simple model that takes into account 
the relativistic tidal force of the central source as well as the Lorentz force 
due to the large-scale electromagnetic field which is assumed to be present 
in the ambient medium. The evolution of the speed of the flow relative 
to the ambient medium is studied. In the force-free case, the resulting 
equation of motion predicts rapid deceleration of the initial flow and an 
asymptotic relative speed with a Lorentz factor of y2. In the presence of 
the Lorentz force, the long-term relative speed of the clump tends to the 
ambient electrical drift speed. 
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1 Introduction 

A significant feature associated with quasars and active galactic nuclei is the 
existence of relativistic outflows known as "jets." The jets frequently come in 
pairs that originate from a central "engine" that is conjectured to be a massive 
rotating black hole surrounded by an accretion disk. Rather similar phenomena 
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are also observed in most of the X-ray binary systems in our galaxy; therefore, 
these are usually called "microquasars." It is thought that the jets are emitted 
along the rotation axis of the Kerr black hole. The double-jet structure is con- 
sistent with the symmetry of a Kerr source under reflection about its equatorial 
plane as illustrated in Figure ^ The observational aspects of the motion of 
astrophysical jets are discussed in QJ [2| E9 Q] . 

The purpose of this paper is to discuss a relativistic tidal effect that can be 
used to predict the speeds of astrophysical jets relative to their ambient media. 
The theory of the formation and propagation of jets per se is beyond the scope 
of this paper. We will consider a Kerr black hole with plasma outflow along its 
rotation axis and study the dynamics of a clump in the flow after it has emerged 
from its source environment. To explore the transient role of the novel relativis- 
tic tidal effect near the source, we assume that large-scale electromagnetic fields 
are dynamically important. We thus ignore the complicated plasma processes 
that actually determine the electrical drift velocity of the ambient medium. Our 
approach may be contrasted with the purely hydrodynamic models of relativis- 
tic outflows: all force-free as well as MHD wind models assume at the outset 
that the plasma flows at the drift velocity. It is important to emphasize the 
idealized nature of the model employed here; our purpose is simply to illustrate 
the regime in which previously ignored relativistic tidal forces are significant for 
the dynamics of relativistic outflows. 

The mechanism for clump formation and the possible values of the clump's 
initial speed as well as other characteristics are beyond the scope of our paper. 
Once emitted, the motion of the clump relative to the ambient medium is subject 
to the ubiquitous tidal forces of the central source together with the Lorentz 
force due to the large-scale electromagnetic field that is assumed to be present 
in the ambient medium. In addition to various other simplifying assumptions 
throughout this work, all radiative as well as nonlinear plasma effects [5] |H] are 
also neglected here. 

With respect to a background Fermi frame that is fixed on the ambient 
medium with coordinates (T, X) such that X = (X, Y, Z) and the Z-axis repre- 
sents the axis of rotation of the black hole, the equation of motion of a plasma 
clump along the rotation axis relative to the ambient medium is given by a 
"Newtonian" equation of motion, involving the superposition of external accel- 
erations, of the form 

d2z „ 

— =A T + A L , (1) 

where At is the relativistic tidal acceleration due to the central rotating mass 
and Al is the Lorentz acceleration due to the electromagnetic field of the am- 
bient medium. The nature of these accelerations are discussed in turn in Sec- 
tions [21 01 and The main dynamical features of the resulting equation of 
motion are described in Section [5J Section contains a brief discussion of our 
results. We employ units such that c = 1 . 
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Figure 1: Schematic diagram representing jet clumps moving relative to the 
ambient medium along the rotation axis of a central black hole. 
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2 The Generalized Jacobi Equation 



Imagine a background global inertial system with coordinates (T, X) and two 
neighboring test particles in free fall in a Newtonian gravitational potential 
$(X). Let Xi(T) and X 2 (T) denote the paths of the two particles. Neglecting 
the gravitational attraction between the test particles, their relative position 
X(T) = Xi(T) - X 2 (T) is determined by the tidal equation 

l^- + Kij (T)X^ = 0, (2) 
where K — [Kij) is a symmetric matrix given by 

Kl ° = dX*dXi (3) 

evaluated along the orbit of one of the particles designated as the reference 
trajectory. Note that tr(i^) = V 2 $ = 4irGpN, where pjsi is the Newtonian 
mass density of the source. If the motion occurs exterior to the source, then the 
matrix K given by J2J is trace-free and harmonic. The tidal acceleration in J2J 
has been evaluated to linear order in the separation of the particles. 

In the theoretical framework of general relativity, the free neighboring test 
particles follow timelike geodesies of the spacetime manifold. Their relative 
motion is described by the geodesic deviation equation. This equation may be 
expressed in terms of a Fermi normal coordinate system that is constructed along 
the reference geodesic. A Fermi frame is an almost inertial coordinate system 
in a finite cylindrical region along the worldline of the reference geodesic \J\. 
We now let (T, X) refer to the Fermi coordinates in such a system, where the 
spatial origin X = is occupied by the fiducial test particle. Then, the Jacobi 
equation (J2J describes the relative motion of the other test particle in Fermi 
coordinates and 

K ij = F R W0 3, (4) 

provided that the speed of relative motion is negligibly small compared to c = 1 . 
Here 

so that the Riemann tensor in the Fermi frame evaluated along the reference 
trajectory (T, 0) is in effect the Riemann tensor projected onto the nonrotating 
orthonormal tetrad A^ that is carried along the fiducial geodesic and upon 
which the Fermi system is constructed. Specifically, the temporal axis of the 
fiducial observer A^ = dx^/dr is the vector tangent to its worldline, r is its 

proper time and A^,i = 1,2,3, are its spatial axes. One can show explicitly 
that in the Newtonian approximation of general relativity (0J reduces to Q , so 
that the formal analogy developed here has a deep physical basis. 

The deviation between geodesies is taken into account only to linear order 
in the Jacobi equation; in fact, higher-order terms can be neglected so long as 
the deviation is very small compared to the radius of curvature of spacetime 1Z. 
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On the other hand, in certain circumstances, the relative speed of neighboring 
geodesies may not be small compared to unity. The corresponding generalization 
of the Jacobi equation is given by 



d 2 X l 
dT 2 



F RmajX : > + 2 F R i kjoV k X : > 



+ (2 F R 0kj0 V i V k + 2 - F R ikji V k V i + 1 F R Qkjt V i V k V t )Xi = 0. (6) 

This generalized Jacobi equation has been discussed in detail in [5]. If we 
limit our attention to motion along a fixed direction in space, say along the 
Z-direction that characterizes the rotation axis of a Kerr system in the Fermi 
frame, then © reduces to 

^- + K (1-2V 2 )Z = Q, (7) 

where V = dZ/dT and k = F Roso3. This equation has an interesting feature 
that was first pointed out in |SJ : Special solutions of |J7J exist that correspond to 
uniform rectilinear motion with limiting velocities of Vt — ±(2)~^- For k < 0, 
these correspond to attractors if 

k(T){T + C a )dT = -co (8) 

To 

for an arbitrary real number Co and an arbitrary initial time Tq [Sj. 

It is useful to consider the application of J7J) to motion along the rotation 
axis of a Kerr source. The main results are expected to hold qualitatively for any 
central rotating and axisymmetric configuration such as, for instance, the Kerr 
spacetime endowed with an infinite set of higher moments [5] . The motion of free 
test particles along the symmetry axis of an exterior Kerr spacetime representing 
the stationary gravitational field of a source of mass M and angular momentum 
Ma is given by 

dr\ 2 2 1 IGMr 

in Boyer-Lindquist coordinates. Here 7 > is an integration constant; in fact, 
in the 7 > 1 case it has the interpretation of the particle's Lorentz factor at 
infinity (r — > 00). We take the reference geodesic to represent the motion of the 
ambient medium with dr/dr > along the Z-axis in the exterior Kerr spacetime 
and we construct a Fermi coordinate system (T, X, Y, Z) in the neighborhood of 
this geodesic such that as the reference geodesic is approached, (T, X) — > (t, 0). 
Evaluating the Ricmann curvature tensor in this frame, we find that n — > k is 
given by 

GMr(r 2 -3a 2 ) 

k - 1 (r 2 + a 2 )3 ■ y™) 

That (|10|) is explicitly independent of 7 is a consequence of the fact that Kerr 
spacetime is of type D in the Petrov classification and hence the axis of rotation 
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provides two special tidal directions for ingoing and outgoing trajectories ^3 
lllj . Moreover, a = corresponds to radial motion in the exterior Schwarzschild 
spacetime. For a Kerr black hole < a < GM, where a = GM corresponds to an 
extreme Kerr black hole. To emphasize that our treatment involves motion along 
the rotation axis of a Kerr system, it is useful to introduce Cartesian coordinates 
(x,y,z) defined in terms of Boyer-Lindquist coordinates (r,6,tp) as usual: x = 
r sin# cos ip, y — r sin 6* sin ^ and z = rcosO. In Section [3] we will employ this 
z-coordinate along the Kerr rotation axis that in the Fermi system becomes the 
Z-axis. Moreover, we will assume that the reference particle at Z = is such 
that 7 = 1 and r > ^/3a along its geodesic path. The 7=1 condition implies 
that the reference particle moves away from the black hole forever, coming to 
rest at infinity. The collection of such particles provides an ambient medium 
with respect to which the motion of the jet clump can be described. That is, 
this ambient medium provides a dynamical way of characterizing the rest frame 
of the black hole. 

Imagine two neighboring test particles at r and r + 5r along the rotation 
axis of a central source. Far away from the source, the Newtonian gravitational 
acceleration of the reference particle at r is —GM/r 2 and for the other particle 
at r + Sr is —GM/(r + Sr) 2 . The relative or tidal acceleration is then given 
by 2GMSr/r 3 , which by {TO} is -kSr for r » GM and r >> a. Therefore, 
(7} and HI U| l reduce to the familiar Newtonian result once the relativistic tidal 
acceleration proportional to V 2 << 1 is neglected. 

The speed of a clump can be measured by observing the displacement of 
the clump in the flow relative to the ambient medium ^|2J03E]. Therefore, 
consider a clump along the axis of rotation moving rapidly past the reference 
particle belonging to the ambient medium. The equation of relative motion is 
given by Q with 

A T = -k(l - 2V 2 )Z, (11) 

where k is determined via equations © and (|10fl . Under our assumption we 
have that r 2 > 3a 2 , and thus k < 0. Also, we have assumed that 7 = 1 
in jnjl for the ambient medium surrounding the black hole; therefore, condition 
||HJl is satisfied, the two special solutions of l|ll|) are indeed attractors and the 
clump speed (predicted by our model) tends to w 0.7 over time In |S] a 
similar assertion was also made for 7 < 1, which turns out to be erroneous in 
general; in any case, it is not relevant for a black hole source with a < GM. 
The qualitative dependence of the solution of the generalized Jacobi equation 
upon 7 > 1 is depicted in Figure [21 m constructing this figure, we consider 
equation (|11|) simply as a nonlinear differential equation and ignore its possible 
physical limitations. 

It is important to recognize two salient features of equation ljll|l for the 
force-free motion of an ultrarelativistic clump relative to the ambient medium 
with 7=1. First, the initial deceleration of the clump can be quite significant 
within the range of validity of the generalized Jacobi equation (|11|) . For instance, 
Figure depicts the deceleration of clumps with initial Lorentz factors of Tq = 
20, 10 and 5, where we have introduced the Lorentz factor of the clump T — 
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Figure 2: Schematic representation of direction fields and trajectories for the 
force-free case given by the generalized Jacobi equation. Equation I jllj l is in- 
tegrated with the initial conditions that at T = 0, r = ro > Via, Z = and 
V = Vq > 0; thus, initially dV/dZ = as well. The top panel depicts a trajec- 
tory for 7 = 1. In this case the velocity is asymptotic to l/v2- The bottom 
panel depicts a trajectory for 7 > 1. In this case, the trajectory is asymptotic to 
some constant velocity, depicted by the dashed line, that is closer to l/y/2 than 
the initial velocity. The limiting velocity can be calculated, in principle, from 
the integration of the generalized Jacobi equation VdV/dZ = —k(l — 2V 2 )Z, 
namely, 1 - 2V 2 = (1 - 2V 2 ) exp(4 k(T(Z'))Z' dZ'). For 7 = 1, the integral 
as Z — > 00 is — 00 in agreement with (jSJ) and hence the limiting velocity is given 
by V = 1/V2. 
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Figure 3: Plot of the Lorentz factor T = (1 - V 2 )~ 1/2 versus T/(GM) based 
on the integration of equation with initial data Tq/{GM) = 10, Z(0) = 0, 
7(0) = 2V6/5 re 0.980, sVU/10 « 0.991 and \/399/20 ps 0.999 corresponding 
to To = 5, 10 and 20, respectively. Here we assume that a/(GM) = 1 and 
7=1. At T/(GM) — 10 the corresponding Lorentz factors are T ps 2.945, 
3.291 and 3.398, respectively. For the sake of comparison, the horizontal line 
indicates the Lorentz factor y/2 corresponding to the terminal speed. Further 
away from the central source, the tidal forces decrease and hence the initial 
deceleration is correspondingly weaker. For instance, with ro/(GM) = 100, the 
clump decelerates from ro = 5 to T = 3 over a time interval given approximately 
by T/(GM) = 300. 
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(1 — F 2 ) -1 / 2 , and r corresponds to the initial clump speed Vq at r — r and 
T = Z = 0. Launched from ro = 10 GM, the clumps decelerate to T ps 3 in 
about T = 10 GM ps 5 x 10 -5 (M/Mq) sec. Although the distance of such a 
clump relative to the reference particle increases rapidly and soon the validity 
of the generalized Jacobi equation breaks down, the subsequent motion can 
be measured with the same model equation relative to the motion of another 
suitable reference particle in the ambient medium. This observation uses the 
second remarkable feature of eauation (|ll|l : it does not explicitly depend on the 
choice of the reference particle in the ambient medium with 7 = 1. To illustrate 
this point, imagine a sequence of n reference particles at r,, i = 1, 2, • • • , n, all 
with 7 = 1 along the rotation axis of the central source such that r\ < T2 < ■ ■ ■ < 
ri < rj+i < ■ ■ ■ < r n . We can set up distinct Fermi frames along the worldlinc 
of each of these reference particles such that within each interval [r^r^+i] the 
generalized Jacobi equation is valid in the Fermi frame based on r^. Once the 
clump reaches n+i, we can switch to the new Fermi frame based on rj+i. So 
long as the initial speed of the clump relative to n+% is above 1/V2, the model 
predicts that the clump continues to decelerate toward the terminal speed. In 
this sense, the speed of the clump relative to the ambient medium approaches 
the terminal speed over time. 

In the absence of any nongravitational forces, equations Q and (fTTf) express 
the geodesic equation for the center of mass of the clump in Fermi coordinates. 
Once the electromagnetic forces are included, the general equation of motion 
takes the form 

, ,d 2 x p „„ dx v dx p , 

1- r£ 

ds 2 p ds ds 

where p' m is the invariant mass density of the clump, F^ v is the Faraday ten- 
sor of the exterior electromagnetic field in the ambient medium and j M is the 
charge current of the clump. Expressed in the Fermi coordinate system (T, X), 
equation (|T2")l reduces to equation QJ, where Al is due to the Lorentz force. In 
evaluating the Lorentz acceleration Al in the following section, we employ the 
physically reasonable approximation that the effects of spacetime curvature can 
be neglected for the sake of simplicity; thus, we treat the Fermi frame established 
along the reference geodesic as an inertial frame of reference in Section |2| 

It is conceivable that the electromagnetic field configuration is such that 
F p a 2° — in equation i|12|) and hence the motion of the clump is force- free (see, 
for example, ^2])- This may be the case, for instance, beyond the launching 
point of the flow, a few gravitational radii from the central source. In the 
force-free case, the detailed analysis of |E] is applicable in the sense explained 
above and the clump speed tends to 1/V2 ~ 0.7, corresponding to a Lorentz 
factor of \/2 « 1.4. Our results may be of interest in connection with Galactic 
superluminal sources 00]; in fact, as already pointed out in some of these 
superluminal jets have speeds that may be near 1/V2 ps 0.7. In the rest of this 
paper, we investigate clump motion that is not force-free. 



P'm{—+K— — )=F\f, (12) 
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3 Lorentz Acceleration 



Imagine a definite clump of plasma with constant mass m moving with respect to 
the ambient medium with velocity V = Vri, where n is a fixed direction in space 
that will be later identified with the rotation axis of the central gravitational 
source. The Lorentz force law for the motion of the clump can be written as 

TTlV 

(6V)(pE + J x B) • n, (13) 



where 8V, p and J are, respectively, the volume, charge density and charge 
current of the clump with respect to the inertial frame (T, X) . We assume that 
in the rest frame of the clump m = p' m 8V', where p' m and 6V are the proper 
density and volume of the plasma clump, respectively, and are assumed to be 
constants throughout the motion under consideration here. Primed quantities 
refer to the rest frame of the clump. It follows from Lorentz contraction that l|13fl 
may be written as 



d ( V \ VI - V 2 



■ 



dT\^i-V 2 
Moreover, we note that for one-dimensional motion 

d / V \ „ ni _, M dV 



(/>E + JxB)-n. (14) 



) = (i-vT 3/2 %, (15) 



dT\^l-V 2 J dT 
so that equation (|14f> reduces to 

^ = (lzp! ( , E + JxB ).a tie, 

dT p' m 

We now turn to the electromagnetic aspects of (|16|l and assume that in the 
rest frame of the clump p' = and that J' = crE' in accordance with Ohm's 
law. The electrical conductivity a is a tensor in general; however, we take a to 
be a scalar for the sake of simplicity. It follows from 

p' = T(p-V-3) (17) 

that p = V • J = VJ\\. This fact immediately implies that J!, = T^ 1 J|| using 
the transformation 

jjj - r(j„ - v P ). (is) 

Thus, we have 

J±=J±, ^||= r ~ lj H ( 19 ) 

and at the same time 

= T(E + V x B) x , E '\\= E W (20) 
It follows from (gHJ and J' = crE' that 

J = crr(E + V x B). (21) 
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Thus p = VJ\\ = <tTVE\\ and 

(pE + J x B) • n = aTVE 2 + aT[(E + J x B) x B]|. (22) 
Using the identity 

[(VxB)xB] ll = -V(B 2 -B 2 ) = -VB 2 ± , (23) 

we find that 

(pE + J x B) ■ n = TaB 2 (w — ^ 11 V), (24) 

where 

(ExB)-n 

w = J (25) 

and 

E x B 

w ^ = -p- ( 26 ) 

is the electrical drift velocity ; 13. . To proceed further, we need a model for the 
electromagnetic field configuration in the medium surrounding the black hole. 
This is considered in the next section. 



4 A Simple Model 

Let us now consider the magnetically dominated environment in which the 
clump is moving. In some models for such an environment |14| and in cylindri- 
cal coordinates (p, cf>, z) with n = z, integrating Maxwell's equation V • B = 
implies that B z ~ p~ 2 . If angular momentum is approximately conserved in the 
electromagnetic field, then B^ ~ p^ 1 . Thus, at large distances from the source, 
jets and collimated winds are often considered to have predominantly toroidal 
magnetic fields. Assuming that the plasma is efficient at maintaining a very 
small proper electric field through its conductive properties, this implies that 
E is approximately radial in cylindrical coordinates |14j . It is clear from this 
brief description that in any reasonably realistic scenario the electromagnetic 
field configuration would be rather complicated. However, we seek a simple 
model situation that would render the resulting equation of motion amenable 
to mathematical analysis. Therefore, we assume that the average electric field is 
primarily radial, E = Ep, and the average magnetic field is primarily azimuthal, 
B = B<fi, such that E < B. In view of the symmetry of the central source about 
its equatorial plane, in this paper we concentrate — for the sake of simplicity — 
only on the clump moving along the positive z-axis. It follows from (|25|l that 
the electrical drift speed is given by w = E/B < 1. Moreover, l|24|) takes the 
simple form T(aB 2 )(w — V). 

It turns out that in the rest frame of the clump E' = E'p and B' = B'<fi, 
i.e. the fields have the same configuration as in the global inertial frame (T, X). 
Thus following the analysis in §2.10 of the conductivity a should be iden- 
tified with a± as the current J' = cr^E' will be crossing the magnetic lines of 
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force in the local rest frame of the clump. On the other hand, for motion along 
the magnetic lines of force, the electrical conductivity un would have its usual 
value 

n P e 2 T r 

°\\ = - J — 1 , (27) 

Me 

where — e, n e and p e are, respectively, the charge, number density and mass 
of the electron in the electron-proton plasma in the rest frame of the clump. 
Here the relaxation time t c = v^ 1 is the average time interval between electron 
collisions and hence v c is the average frequency of electron collisions per second. 
It follows from the discussion in §2.10 of that in case of a tenuous plasma 
in a strong magnetic field such that i/j < fl', we have 



where Q e = eB' / fi e is the cyclotron frequency for a free electron. Thus 

<r ± B 2 fl i e \fB\* 



where fx p is the proton mass and p' m = n p n p + n e p e = n e (fi p + fi e ) ~ n e ^, p , since 
p' = implies that the density of protons in the clump is the same as that of 
electrons, n p = n e , and [i p j '/x e w 1836 3> 1. Moreover, 

Bj_ = T(B — V X E)j_, (30) 

which implies that 

^ = r(i-wV). (3i) 

Combining these results with l|16|l reveals that the Lorentz acceleration is given 

by 

dV (1 ~ V 2 f/ 2 

J ' = s=°Vw ( "" n (32) 



(t^W. (33) 



where 

a = 

\p P > 

It is interesting to note that for copper at room temperature v c ~ 10 14 sec _1 . 

It follows from equations and (|33|l that the Lorentz acceleration Al 
vanishes if the clump speed is equal to the drift speed. Moreover, Al is propor- 
tional to the collision frequency v c . This is due to the fact that only through 
collisions can an appropriate interior current be established that could lead to 
the acceleration of the clump along the flow direction. 

It remains to provide an estimate of v c for the case of the electron-proton 
plasma under consideration here. This is a rather complicated problem and 
wc therefore limit our investigation to the simpler case of collisions among the 
electrons while neglecting the motion of the protons. Thus we assume that r c 
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can be estimated using the "self-collision time" of electrons given by equation 

(5-26) of ua 

, 0.266 r 3 / 2 

t c = — — sec, (34) 

n e in A 

where T is the absolute temperature and n e is the number of electrons in the 
clump per cm 3 . Here A is essentially the ratio of two lengths: the Debye shield- 
ing distance over the characteristic impact parameter of electron collisions such 
that the deflection angle in the orbital plane is equal to 7r/2. One can use Table 
5.1 of to find the values of In A in terms of T and n e . According to this 
table, for a fixed T, In A decreases very slowly but monotonically as n e increases 
from n e = 1 cm~ 3 to n e = 10 24 cm~ 3 ; for example, for T = 10 8 K, In A mono- 
tonically decreases from 34.3 for n e = 1 cm~ 3 to 6.69 for n e = 10 24 cm~ 3 . Thus 
for T = 10 8 K and n e = 10 6 cm -3 , we find In A = 27.4 and hence from (J31J with 
v c (t*) -1 we obtain a » 0.56 x 10" 7 sec~ 1 . For fixed T, a increases almost 
linearly with n e such that for n e = 10 15 cm -3 in the case under discussion, we 
have a « 35 sec -1 [T3] . These considerations illustrate the fact that a could 
have a considerable range of values depending upon the electron temperature T 
and density n e . 



5 Equation of Motion 

The equation of motion of the clump is given by 

% + fc ™ - ^ Z - ai ir^w (w - n (35) 

where V — dZ/dT, a is a characteristic of the plasma clump, w is a characteristic 
of the electromagnetic field of the ambient medium and k is essentially the 
curvature of the central source given by 

2GMz{z 2 - 3a 2 ) 
K T ) = , 2 , 2 x 3 , 36) 

while z(T) is determined using the geodesic equation for the background 

dz 2 2 2GMz 

fe) =7 - 1 + ^T^- (37) 



The equation of motion (|35(l has been derived for the clump moving along 
the positive Z-axis. Recall that z in equations l|36l) and 1)3 7|) is a substitute 
for the radial coordinate r. Therefore, let us note that equations l|35l) - (|37|l are 
invariant under the transformations Z — > —Z. w — > — w, V — > —V and z — > z. 
Hence our results hold for the clump moving along the negative Z-direction as 
well. 

It is useful to put equations l|35 |l - H37|l in dimensionless form. To this end, 
we assume that all lengths are measured in units of GM, which is one- half of 
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Figure 4: The top panel depicts the graph of the Lorentz factor T, given by 
r(T) := (1 — y(T) 2 ) -1 / 2 , versus time specified in hours. The bottom panel 
depicts the graph of T versus time specified in seconds. For these graphs, a = 1, 
7 = 1, a = 10~ 7 , w = V3/2 ps .866, r /(GM) = z(0) = 100, Z(0) = and 
V(0) = 2V6/5 re 0.980. We assume that M = 10Af Q , hence time is measured 
here in units of GM = 10 GM R pa 5 x 10~ 5 sec. 



the gravitational radius of the central source. Defining dimensionless quantities 
k, a and a as 

k = (GM) 2 k, a = GMa, a = (38) 
v ' GM K ' 

we recover equations (|35[) - (|37H in dimensionless form once we let k — > k, a — > a, 
a^a and GM — ► 1 . 

The equations of motion in dimensionless form are equivalent to the dynam- 
ical system 

dz ( r, 2z \V2 



dZ_ 

dT 



V. 



dV ^{l-V 2 f/ 2 t . 2z(z 2 -3a 2 ), 1 „ 

where we assume that 7 > 1. 

We have proved the following result: If a > 0, < w < 1, z(0) > V3a, and 
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< V(0) < 1, then liniT^oo V(T) = w. That is, as time increases the clump 
speed approaches the drift speed of the ambient medium. 

Equation (|35|l and the system (|39|l show how the tidal equation can in prin- 
ciple be generalized to include nongravitational forces. On the other hand, it 
is necessary to emphasize the qualitative significance of equation (|35fl , since its 
gravitational part is valid only within a distance 1Z of the fiducial particle and 
its electromagnetic part is based on a rather simple model. If the tidal part is 
completely ignored, then the Lorentz force would lead to a monotonic deceler- 
ation of an initially ultrarelativistic clump toward w. Numerical experiments 
based on the system l|39[) demonstrate that the tidal part initially dominates, 
causing a very rapid drop in T toward -\/2, and this occurs mostly within the 
domain of validity of the tidal term, but then the electromagnetic term takes 
over and T slowly tends toward (1 — w 2 ) -1 / 2 . Figure 0] depicts a profile for 
the graph of T versus time, where the parameter values and initial conditions 
are chosen to be typical for a microquasar with M — IOMq. According to our 
model <|39[l with electrical drift speed corresponding to T = 2 and with Lorentz 
acceleration coefficient a = 10~ 7 , an initial Lorentz factor of T — 5 decreases to 
F pa 1.5 in about 0.25 sec and then very slowly increases to within 1% of F = 2 
after 48 hours. More generally, for a < 10~ 2 the timescale for the initial rapid 
decrease in T is about 0.1 sec. The corresponding total relaxation time to T = 2 
is approximately (50a) _1 sec, so that for a ranging from 10~ 3 to 10 -11 , it ranges 
from 20 seconds to 64 years. For a quasar with M — 10 8 A/ Q , the timescale for 
the initial rapid decrease in T would be about two weeks. 

To see that our result holds in general, let us note that dz/dT > for T > 0. 
Thus, z is an increasing function of T, which is defined for < T < oo. By 
viewing z as a new "time" in the dynamical system . it suffices to consider 
the system 



dZ 1 



■V. 

-2n 



dz R(z) 

« * ( s (izZ^ ( „_ y) + 2£Ml (1 _ 2nz ), (40) 

dz R(z)\ (1-wVy (z 2 + a) 3 > 



where 



R{z) := ( 7 2 - 1 



2z \V2 



z z + a ' 

In this system dZ/dz > 0. Hence, as before, it suffices to consider the scalar 
differential equation 

dV ^ {1-V*ff\ nr7 S(Z)(S(Z)i-3a 2 )(l-2V*) 

dz = a { i- wV? v {w - v) + 2Z {S (zr + a*yv ' (41) 

where S is the function such that S(Z(z)) = z. 

Direction fields of dV/dZ in the (Z, V)-plane are depicted in Figure|S]for the 
two cases w > l/v2 and w < 1/V2. Graphs of typical solutions are also drawn. 
Since V stays bounded and the curvature k(T) approaches zero very fast, it can 
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Figure 5: Schematic representation of direction fields and trajectories for the 
dynamical system l(Hjl. 



be shown that the second term on the right-hand side of equation (|41|) rapidly 
approaches zero as T — > oo. Also, with this term set to zero, the dynamical 
system (|41|) has an asymptotically stable steady state at V = W. For large Z, 
the first term is dominant; therefore, the limiting value of the solution is w. 

6 Discussion 

This paper is devoted to a simple model of the evolution of the speed of a 
clump in a relativistic flow once it emerges from the environment immediately 
surrounding the central source. Combined with a realistic treatment of plasma 
effects, our theoretical approach should be of interest in the study of astro- 
physical jets. If the motion of the clump is force-free, then it follows from the 
repeated application of the generalized Jacobi equation that the clump speed 
tends to Vr = 1/V% corresponding to T = \[2. On the other hand, if the force- 



16 



free condition is not satisfied, the clump speed tends to Vl = w, where w is 
the electrical drift speed. In fact, if the tidal component initially dominates, as 
would be the case for most microquasars, then starting with an initial speed al- 
most equal to unity, the tidal force is responsible for the initial rapid decrease of 
the velocity toward 1 / \/2, but then the Lorentz term takes over and the velocity 
slowly approaches w. On the other hand, if the Lorentz term is dominant, as 
would be the case for most quasars, there is in effect a relatively slow decrease 
of the velocity toward w. 
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